Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1155 1791 (N56E02 L= 93)  1387    287        0 (1034/1760-1298/1822)  -1.9 E16-W14 N54-N57
      2  480 1605 (N30E52 L= 42)  2281    369        3 ( 383/1521- 569/1670)  21.8 E60-E43 N24-N35
      3 1217 1572 (N38W07 L=101)   818    188        0 (1129/1554-1303/1593)   8.7 W00-W13 N36-N40
      4  283 1519 (N22E70 L= 24)  1201    108        3 ( 249/1473- 307/1565) -56.3 E76-E66 N18-N24
      5  385 1452 (N19E53 L= 42)   631    104        1 ( 364/1403- 407/1501) -60.3 E55-E48 N15-N22
      6 1752 1343 (N29W52 L=147)  3814    409        9 (1549/1279-1911/1414)  -3.7 W33-W80 N25-N31
      7  256 1327 (N10E61 L= 34)   418     86        0 ( 244/1281- 273/1368) -80.5 E62-E59 N06-N12
      8  187 1212 (N02E66 L= 29)  1697    179        1 ( 142/1131- 215/1269) -31.0 E75-E61 S02-N05
      9  779 1045 (S01E15 L= 80)   645    225        2 ( 756/ 943- 805/1148) -85.2 E16-E13 S07-N05
     10 1890 1106 (N15W67 L=162)   503     79        2 (1873/1056-1910/1133)  68.2 W64-W70 N12-N17
     11  473  998 (S07E35 L= 60)   880    271        1 ( 412/ 895- 537/1119)  74.1 E39-E31 S14-N00
     12 1345 1094 (N09W19 L=114)   402    136        1 (1293/1065-1396/1118)   8.1 W15-W22 N06-N10
     13 1073  989 (S01W03 L= 98)   480    115        1 (1025/ 945-1108/1031)  53.1 W00-W05 S03-N01
     14 1341  658 (S18W25 L=120)   427    166        0 (1268/ 631-1398/ 703)  26.6 W20-W28 S20-S14
     15 1396  549 (S24W32 L=127)   399     89        2 (1349/ 529-1444/ 573)  31.0 W28-W34 S25-S22
     16 1131  387 (S39W18 L=113)   657     75        2 (1089/ 367-1175/ 415)  38.7 W15-W21 S41-S36
  

    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