Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1090 1573 (N38E00 L=298)   324     90        0 (1051/1556-1136/1594)   0.0 E03-W03 N36-N39
      2 1235 1525 (N36W11 L=310)   267     73        0 (1208/1498-1269/1539)  14.0 W09-W13 N33-N37
      3 1698 1485 (N37W51 L=349)   893    191        2 (1628/1438-1768/1518) -22.6 W44-W57 N33-N38
      4  527 1450 (N24E39 L=259)   315     51        0 ( 506/1421- 552/1479)  76.0 E40-E38 N22-N26
      5  632 1439 (N24E31 L=268)   547    151        0 ( 572/1404- 695/1466) -15.3 E36-E25 N22-N25
      6 1554 1340 (N25W35 L=333)  1928    427        5 (1371/1234-1692/1458) -28.4 W21-W45 N19-N32
      7  228  983 (S08E58 L=241)  1348    216        1 ( 200/ 875- 250/1094) -69.0 E62-E55 S14-S01
      8  736  738 (S18E16 L=282)   291     80        1 ( 704/ 713- 780/ 760)  38.7 E18-E13 S20-S16
      9 1790  662 (S16W62 L=  1)   952     94        2 (1761/ 627-1813/ 708)  90.0 W58-W65 S17-S12
  

    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