Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  905 1428 (N31E06 L=120)   689    233        0 ( 843/1339- 946/1511) -67.4 E10-E03 N24-N36
      2 1291 1343 (N23W20 L=146)   753    218        0 (1211/1290-1365/1373)   0.0 W13-W24 N19-N24
      3  262 1292 (N23E59 L= 68)  1086    138        2 ( 241/1214- 289/1350)  82.9 E61-E55 N18-N26
      4 1835 1275 (N13W65 L=192)   896     76        2 (1799/1241-1868/1308) -20.6 W61-W70 N11-N15
      5  192 1234 (N19E67 L= 60)  1431    123        3 ( 166/1180- 219/1281)  66.8 E71-E63 N15-N22
      6 1648 1181 (N10W43 L=170)  1193    270        6 (1550/1132-1772/1231) -21.3 W36-W54 N05-N13
      7 1858  941 (S07W63 L=189)   522     40        0 (1840/ 915-1876/ 976)  90.0 W60-W65 S08-S05
      8 1539  838 (S10W33 L=159)   416     86        2 (1520/ 795-1564/ 879)  68.2 W31-W34 S12-S07
      9 1590  744 (S17W37 L=163)  1799    316        6 (1494/ 651-1671/ 844)  42.5 W30-W43 S22-S10
     10  974  707 (S15E05 L=121)   671    218        1 ( 928/ 612-1040/ 800)  54.0 E08-W00 S20-S09
     11  520  738 (S11E35 L= 91)   316     75        1 ( 493/ 691- 542/ 762)  68.2 E36-E32 S14-S09
     12 1816  692 (S23W62 L=188)  1314    122        2 (1766/ 656-1867/ 734)  14.0 W56-W69 S25-S20
     13 1599  664 (S22W39 L=165)   276     59        0 (1573/ 651-1631/ 678)  11.3 W36-W41 S23-S21
     14 1540  612 (S25W34 L=161)   407     79        1 (1512/ 580-1567/ 645) -33.7 W31-W37 S27-S22
     15  578  486 (S29E36 L= 90)  2408    374        8 ( 434/ 409- 669/ 602)  24.0 E53-E26 S34-S20
     16  785  514 (S27E20 L=107)   322     77        2 ( 757/ 486- 810/ 556) -68.2 E21-E17 S29-S24
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y