Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  646 1700 (N61E07 L= 58)  2066    454        5 ( 472/1600- 814/1819)   5.6 E31-W20 N56-N64
      2 1669 1359 (N06W50 L=115)   315     56        2 (1642/1337-1689/1384)  14.0 W47-W52 N04-N07
      3 1207 1349 (N19W20 L= 85)   558    158        2 (1125/1312-1284/1366) -29.1 W14-W24 N16-N22
      4 1840 1070 (S15W55 L=121)  3731    446        9 (1789/ 882-1910/1290)  48.0 W46-W64 S24-S02
      5 1898 1167 (S12W66 L=131)  2004    292        6 (1859/1028-1945/1268) -79.4 W61-W71 S21-S05
      6 1493 1183 (N03W31 L= 96)   280     89        0 (1479/1135-1507/1232)  53.1 W29-W33 N00-N04
      7  243 1126 (N30E54 L= 12)   457     99        0 ( 200/1108- 289/1150) -15.3 E59-E48 N29-N32
      8 1746 1025 (S14W45 L=110)   935    145        4 (1713/ 958-1788/1090)  32.0 W40-W49 S17-S10
      9  174 1047 (N27E63 L=  2)   405     82        0 ( 136/1025- 212/1070) -15.3 E69-E57 N24-N28
     10  607 1005 (N16E25 L= 41)   710    347        4 ( 463/ 954- 741/1053)  -9.0 E34-E15 N12-N19
     11  856  973 (N08E11 L= 55)   376    171        1 ( 788/ 940- 906/1024) -60.3 E13-E07 N05-N12
     12  253  832 (N13E57 L=  8)   914    137        4 ( 221/ 755- 278/ 885)  74.1 E60-E53 N08-N15
     13  810  682 (S07E21 L= 44)   289     55        2 ( 777/ 663- 839/ 699) -45.0 E23-E19 S08-S05
     14 1727  588 (S42W39 L=105)  2104    221        4 (1647/ 534-1786/ 674)  14.9 W30-W45 S46-S35
     15  583  647 (S04E36 L= 29)   312     47        2 ( 554/ 626- 616/ 673) -11.3 E38-E33 S05-S02
     16  469  397 (S17E57 L=  8)  2010    295        7 ( 404/ 326- 581/ 459) -51.3 E63-E50 S23-S12
     17 1251  344 (S42E07 L= 59)  1163    175        3 (1177/ 306-1326/ 384)   0.0 E12-W00 S43-S38
     18  945  179 (S45E40 L= 26)  1092    114        3 ( 883/ 155-1002/ 212)  -4.8 E46-E33 S47-S42
  

    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