Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1325 1758 (N56W19 L= 85)   847    129        2 (1256/1736-1394/1777)   0.0 W10-W27 N54-N58
      2  723 1530 (N29E28 L= 39)  1141    256        2 ( 688/1397- 782/1604) -81.9 E30-E24 N20-N34
      3  405 1530 (N25E55 L= 12)   577     97        0 ( 378/1487- 445/1565)  59.0 E56-E52 N21-N27
      4  306 1469 (N20E63 L=  4)   552    105        0 ( 263/1430- 337/1495) -11.3 E71-E56 N17-N21
      5  296 1410 (N16E60 L=  7)   540    126        1 ( 246/1391- 340/1437)   9.5 E66-E54 N14-N17
      6 1715 1150 (N16W46 L=113)   480    118        0 (1675/1101-1739/1205) -51.3 W43-W48 N12-N18
      7  461 1154 (N03E38 L= 29)   738    192        0 ( 389/1085- 509/1201) -11.3 E43-E33 S01-N05
      8 1534  954 (N02W33 L=100)   545     81        1 (1500/ 908-1556/ 993)  90.0 W30-W34 N00-N04
      9  906  912 (S07E06 L= 61)   478    164        0 ( 861/ 848- 960/ 967)  58.0 E08-E03 S11-S03
     10 1742  649 (S15W58 L=125)  1635    333        3 (1621/ 567-1821/ 737)  42.7 W47-W64 S20-S08
     11 1699  496 (S25W63 L=129)   609     82        0 (1667/ 473-1732/ 523)  45.0 W57-W66 S26-S22
     12 1436  318 (S41W48 L=115)  1369    188        1 (1363/ 280-1522/ 379)  41.6 W41-W52 S44-S35
  

    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