Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1324 1613 (N37W01 L=111)  1842    292        1 (1207/1488-1437/1719) -25.5 E10-W12 N31-N41
      2  605 1562 (N13E41 L= 69)   597    143        1 ( 542/1533- 681/1602)  40.6 E44-E37 N10-N17
      3  252 1462 (N01E69 L= 41)  1150    118        2 ( 223/1406- 290/1514)  71.6 E71-E66 S02-N04
      4  919 1475 (N16E18 L= 92)   310     87        0 ( 874/1454- 961/1491)   0.0 E21-E15 N14-N16
      5 1418 1409 (N26W13 L=123)   302     44        1 (1392/1379-1436/1432) -56.3 W11-W14 N24-N27
      6 1005 1291 (N07E08 L=102)   309     74        1 ( 972/1256-1040/1328) -21.8 E10-E05 N05-N08
      7  258 1123 (S19E53 L= 56)   415    109        0 ( 242/1067- 276/1178) -81.9 E55-E51 S22-S15
      8 1328  598 (S22W31 L=141)   615    147        1 (1252/ 581-1400/ 620)  21.8 W25-W35 S23-S19
      9 1519  393 (S26W55 L=164)   737    130        0 (1467/ 373-1583/ 413)  26.6 W50-W59 S28-S22
     10 1641  356 (S23W68 L=178)   922    177        2 (1596/ 331-1692/ 389)  45.0 W63-W72 S25-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