Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1106 1776 (N54E08 L=163)   750    179        2 (1028/1750-1185/1801)  23.2 E15-W00 N50-N58
      2 1343 1749 (N57W18 L=189)   435     53        1 (1321/1719-1371/1771)  23.2 W15-W22 N53-N59
      3 1438 1725 (N56W30 L=201)   589     44        1 (1415/1696-1465/1751)  56.3 W26-W34 N53-N59
      4  315 1493 (N20E64 L=106)  2509    268        5 ( 274/1392- 349/1609) -47.1 E73-E58 N13-N27
      5  721 1559 (N30E30 L=140)   430    115        3 ( 702/1499- 743/1603) -63.4 E32-E28 N25-N33
      6  242 1451 (N16E73 L= 98)   945    107        0 ( 218/1419- 272/1503)  71.6 E75-E69 N13-N19
      7 1407 1442 (N32W21 L=192)   470    136        0 (1351/1407-1466/1467) -11.3 W16-W26 N30-N33
      8 1783 1346 (N31W57 L=228)   550     70        0 (1746/1328-1815/1366)  15.3 W52-W63 N29-N32
      9  872 1322 (N17E14 L=157)   197     66        0 ( 840/1303- 909/1337)   0.0 E15-E10 N15-N17
     10  684  972 (S07E20 L=150)   320     71        1 ( 653/ 938- 709/1005) -45.0 E22-E18 S08-S05
  

    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