Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1095 1615 (N36W21 L=176)  1093    201        0 (1015/1561-1190/1648)  -4.4 W14-W27 N33-N37
      2  614 1575 (N46E16 L=139)   271     82        1 ( 591/1549- 643/1604) -38.7 E18-E13 N44-N48
      3 1754 1371 (N05W57 L=213)   982    132        1 (1702/1317-1802/1414) -54.5 W53-W60 N01-N08
      4  206 1258 (N33E55 L=101)  1424    461        5 ( 121/1208- 283/1318) -25.3 E65-E46 N29-N38
      5 1282  959 (S06W13 L=168)  1254    354        2 (1183/ 779-1347/1122)  41.4 W03-W20 S14-N01
      6 1768 1080 (S11W49 L=204)   309     78        1 (1731/1056-1794/1095) -21.8 W46-W51 S12-S09
      7 1286  707 (S21W08 L=164)   195     56        1 (1254/ 695-1315/ 724)  14.0 W06-W10 S22-S20
      8  901  437 (S29E23 L=133)   803    172        1 ( 876/ 357- 925/ 514)  61.4 E26-E19 S34-S23
      9 1312  456 (S39W05 L=160)   431     92        0 (1268/ 420-1343/ 501)  33.7 W01-W08 S40-S35
     10  519  384 (S24E51 L=104)  1275    340        5 ( 471/ 292- 557/ 491)  65.0 E57-E46 S31-S16
  

    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