Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1433 1683 (N33W48 L= 76)  2037    324        5 (1286/1618-1555/1732) -40.9 W37-W56 N26-N39
      2  506 1672 (N56E33 L=355)  2634    191        9 ( 442/1601- 573/1724)  19.7 E45-E23 N50-N61
      3  776 1556 (N40E05 L= 23)  1244    243        3 ( 686/1486- 886/1600) -47.5 E11-W01 N32-N44
      4  878 1340 (N23E03 L= 25)   471     62        2 ( 849/1309- 918/1372)  14.0 E04-W00 N21-N25
      5  287 1195 (N26E49 L=339)  1451    265        4 ( 237/1100- 339/1320) -61.9 E56-E43 N19-N34
      6 1000 1234 (N14W03 L= 31)   449    141        0 ( 937/1203-1056/1279) -54.5 W00-W05 N11-N18
      7  406 1137 (N20E38 L=350)   836    156        3 ( 380/1066- 440/1205)  69.4 E40-E34 N16-N24
      8 1808 1076 (S11W54 L= 82)  1177    201        0 (1782/ 971-1836/1158)  55.0 W50-W58 S16-S06
      9  159 1014 (N17E65 L=323)  1725    238        3 (  90/ 956- 246/1109) -24.6 E79-E54 N12-N23
     10  311 1044 (N16E47 L=341)  1076    154        4 ( 259/ 992- 374/1085) -33.7 E51-E41 N12-N18
     11 1631  856 (S20W35 L= 63)   380    117        1 (1585/ 821-1689/ 883) -32.0 W31-W40 S23-S17
     12  293  762 (N00E54 L=333)   392     58        0 ( 261/ 740- 323/ 789) -36.9 E57-E51 S01-N02
     13  874  566 (S22E19 L=  9)   365    144        0 ( 832/ 501- 907/ 630)  41.2 E23-E15 S25-S18
     14 1143  566 (S28E02 L= 26)   401    102        2 (1105/ 520-1178/ 608)  45.0 E05-W00 S30-S24
     15  355  570 (S12E56 L=331)   347     57        1 ( 335/ 541- 372/ 598)  53.1 E58-E54 S14-S10
     16  396  479 (S18E58 L=330)  1258    132        5 ( 331/ 442- 455/ 526) -29.1 E63-E53 S21-S14
  

    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