Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  969 1911 (N64W44 L= 36)   480     42        1 ( 948/1882- 984/1935)  18.4 W37-W52 N60-N66
      2  853 1862 (N64W24 L= 16)  1569    199        3 ( 786/1802- 922/1906)  11.3 W12-W38 N59-N67
      3  701 1771 (N62E01 L=352)   445     67        0 ( 682/1739- 726/1803)  35.5 E03-W04 N58-N64
      4  904 1611 (N42W09 L=  2)  1359    313        2 ( 767/1578-1061/1643) -21.8 E01-W19 N35-N45
      5 1675 1538 (N15W59 L= 52)   604     72        0 (1634/1517-1710/1563)   0.0 W54-W64 N14-N16
      6 1186 1496 (N26W23 L= 15)   486     76        1 (1135/1465-1220/1531)   9.5 W19-W25 N24-N28
      7 1829 1418 (N04W69 L= 61)  1389    101        5 (1801/1365-1865/1476)  74.1 W64-W74 N00-N07
      8  365 1411 (N43E40 L=312)   476     56        0 ( 339/1387- 400/1431)   0.0 E43-E36 N41-N44
      9 1491 1299 (N08W34 L= 26)   598    141        0 (1451/1234-1524/1369) -90.0 W32-W35 N03-N12
     10 1821 1282 (S02W60 L= 52)   401     74        1 (1801/1248-1848/1306)   9.5 W57-W64 S04-S01
     11 1154 1082 (N04W08 L=  1)   425    129        0 (1113/1016-1193/1133) -90.0 W06-W09 N00-N07
     12  146 1022 (N23E66 L=287)  5545    255       12 (  82/ 898- 227/1130) -78.7 E79-E55 N14-N29
     13  263 1057 (N23E51 L=302)   367     50        2 ( 237/1038- 292/1078) -31.0 E53-E47 N20-N23
     14  388  908 (N12E42 L=310)   793    192        2 ( 318/ 861- 463/ 977) -45.0 E46-E37 N07-N16
     15 1157  785 (S12W02 L=354)   462    166        0 (1131/ 718-1184/ 876)  53.1 E01-W05 S15-S07
     16 1737  584 (S42W43 L= 35)   639     86        0 (1692/ 566-1784/ 609) -24.8 W37-W50 S45-S38
     17  542  514 (S14E44 L=309)   629    129        0 ( 515/ 451- 573/ 581)  37.9 E48-E38 S17-S10
     18 1558  501 (S42W24 L= 16)   593    135        0 (1502/ 458-1610/ 540)  16.7 W18-W28 S43-S40
  

    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