Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  361 1353 (N35E48 L= 64)   442     61        0 ( 326/1337- 395/1367)  -7.1 E52-E44 N33-N35
      2 1304 1295 (N18W22 L=134)   388    139        0 (1249/1242-1346/1334) -60.3 W19-W23 N14-N21
      3 1231 1251 (N17W16 L=128)   435    202        0 (1197/1179-1291/1331) -78.7 W14-W18 N11-N21
      4  306 1052 (N16E50 L= 62)   324     60        2 ( 289/1022- 332/1086)  56.3 E52-E47 N14-N17
      5 1711 1039 (S04W46 L=158)   475    103        0 (1668/1013-1758/1065) -29.7 W42-W50 S06-S02
      6  985  638 (S17E08 L=104)   578    182        3 ( 930/ 568-1030/ 715)  48.8 E12-E04 S20-S12
      7 1191  497 (S30W03 L=115)   665    236        0 (1107/ 414-1274/ 564)  31.6 E04-W09 S34-S26
      8  640  487 (S23E35 L= 77)   720    197        0 ( 566/ 432- 709/ 552)  25.0 E43-E28 S25-S18
      9  766  408 (S30E29 L= 83)  1033    268        1 ( 661/ 346- 867/ 490)  20.9 E40-E19 S33-S25
     10 1234  406 (S38W05 L=117)   461    111        1 (1182/ 374-1284/ 440)  18.4 W00-W09 S39-S35
  

    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