Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  959 1513 (N38W08 L=126)  1255    308        3 ( 925/1359-1004/1633)  48.8 W01-W15 N28-N45
      2 1185 1600 (N37W27 L=145)   398     57        0 (1157/1570-1219/1622)  18.4 W24-W30 N35-N38
      3 1404 1578 (N30W42 L=159)   595    104        1 (1361/1548-1463/1600)   0.0 W37-W47 N28-N31
      4 1648 1455 (N16W54 L=171)   391     84        1 (1615/1432-1687/1477) -26.6 W50-W57 N14-N17
      5 1018 1283 (N22W05 L=123)   678    145        2 ( 969/1230-1070/1340) -71.6 W02-W07 N17-N26
      6  156 1138 (N28E69 L= 49)  2377    163        6 ( 114/1069- 189/1231)  84.3 E77-E63 N23-N33
      7  893  988 (N08E08 L=109)   411    184        3 ( 819/ 958- 968/1013) -26.6 E12-E04 N04-N09
      8  237  885 (N12E58 L= 59)  1263    283        1 ( 188/ 817- 281/ 966)  53.1 E64-E52 N08-N16
      9 1774  815 (S24W47 L=165)  2321    320        2 (1687/ 671-1836/ 949)  42.9 W38-W54 S31-S16
     10 1345  816 (S12W14 L=132)   372    120        1 (1306/ 774-1383/ 876)  26.6 W11-W17 S14-S09
     11  668  771 (S01E26 L= 91)   286    111        0 ( 655/ 716- 678/ 830)  59.0 E28-E24 S03-N02
     12 1339  715 (S18W12 L=130)   688    182        2 (1290/ 643-1409/ 760)  16.7 W07-W17 S21-S14
     13 1471  684 (S23W20 L=138)   758    301        2 (1387/ 576-1554/ 756)  28.3 W13-W27 S28-S19
     14 1028  620 (S17E09 L=109)   340    172        0 (1004/ 563-1071/ 682)  83.7 E09-E05 S21-S12
     15 1204  659 (S18W03 L=120)   294     50        1 (1170/ 635-1238/ 679) -45.0 W00-W04 S19-S16
     16  876  585 (S16E19 L= 99)   211     37        0 ( 855/ 562- 896/ 612) -90.0 E19-E17 S17-S13
     17  470  445 (S18E53 L= 65)  1342    211        1 ( 375/ 409- 548/ 503)   5.4 E64-E43 S20-S14
     18 1044  326 (S38E17 L=101)  4490    688        5 ( 784/ 211-1277/ 451)   2.7 E42-W02 S42-S32
     19 1667  405 (S53W44 L=162)   581     41        1 (1642/ 389-1696/ 422) -21.0 W38-W51 S55-S50
  

    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