Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  895 1689 (N50W16 L= 97)   408     96        0 ( 869/1650- 942/1718)  12.5 W12-W21 N47-N51
      2 1658 1627 (N19W69 L=150)  1812    165        3 (1612/1546-1701/1683)  68.2 W64-W75 N13-N23
      3 1033 1422 (N28W12 L= 93)   319    127        0 ( 975/1386-1081/1455) -63.4 W09-W14 N25-N31
      4  476 1360 (N40E29 L= 52)   634    164        0 ( 450/1293- 511/1442) -74.7 E31-E26 N35-N46
      5  788 1264 (N26E08 L= 73)   524    161        0 ( 702/1244- 861/1284) -33.7 E12-E03 N23-N29
      6  407 1192 (N31E37 L= 44)   290     77        1 ( 381/1161- 423/1225) -90.0 E39-E35 N29-N33
      7 1874 1150 (S13W63 L=144)  1675    174        1 (1804/1091-1948/1191)  -2.0 W52-W81 S15-S10
      8  450 1137 (N27E34 L= 47)   461    104        3 ( 399/1109- 499/1168)  -6.3 E38-E29 N25-N28
      9  987 1114 (N12W00 L= 81)   200     62        0 ( 955/1090-1014/1141) -76.0 W00-W01 N10-N14
     10  541 1075 (N21E28 L= 53)   347     77        1 ( 523/1043- 566/1108)  56.3 E29-E26 N19-N23
     11  310  982 (N21E49 L= 32)   437     83        0 ( 272/ 945- 353/1014) -39.8 E52-E45 N17-N23
     12  612  938 (N12E26 L= 55)   255     96        0 ( 582/ 900- 649/ 975) -63.4 E27-E24 N08-N14
     13  490  876 (N11E36 L= 45)   698    105        3 ( 463/ 826- 528/ 926)  80.5 E38-E33 N08-N14
     14 1528  761 (S23W23 L=104)  1045    393        0 (1385/ 655-1644/ 901)  14.6 W11-W34 S25-S18
     15 1368  802 (S15W14 L= 94)   293     44        2 (1341/ 776-1396/ 821)   0.0 W12-W15 S16-S13
     16 1788  752 (S33W46 L=127)   692     99        1 (1744/ 715-1826/ 783)   5.7 W41-W51 S34-S31
     17 1220  618 (S22W00 L= 81)   405    183        0 (1163/ 564-1278/ 696)  29.1 E04-W05 S24-S18
     18  714  367 (S24E41 L= 40)  1447    343        2 ( 615/ 249- 801/ 460)  16.5 E58-E31 S31-S19
  

    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