Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1261 1601 (N41W07 L=158)   679    155        1 (1209/1549-1317/1639)  45.0 W03-W12 N36-N44
      2  548 1410 (N16E38 L=113)   484    137        4 ( 497/1383- 589/1446)  59.0 E40-E34 N13-N18
      3  728 1283 (N11E22 L=128)   909    237        5 ( 673/1218- 789/1341)  66.0 E25-E19 N06-N15
      4  143 1065 (S11E69 L= 81)  2800    332        5 ( 114/ 965- 171/1189)  90.0 E73-E66 S17-S03
      5  431 1015 (S09E38 L=112)   247     95        0 ( 389/ 993- 468/1045)  26.6 E41-E35 S10-S07
      6  743  917 (S11E16 L=135)   789    217        3 ( 711/ 827- 785/1010) -59.0 E17-E11 S15-S05
      7 1236  992 (N01W13 L=164)   247     47        0 (1205/ 969-1265/1010) -14.0 W11-W15 N00-N01
      8 1429  913 (S01W27 L=177)   373     96        1 (1389/ 888-1482/ 936)  31.0 W24-W29 S02-N01
      9  292  594 (S40E53 L= 98)   705    157        2 ( 263/ 567- 325/ 621)  26.6 E56-E48 S41-S36
     10 1305  428 (S33W32 L=182)   262     49        0 (1278/ 407-1335/ 446)  56.3 W30-W33 S35-S31
     11 1410  377 (S35W42 L=193)   553    131        1 (1373/ 354-1458/ 409) -12.5 W38-W47 S36-S33
  

    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