Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  571 1676 (N20E51 L=132)  1728    309        6 ( 516/1613- 639/1741) -24.8 E59-E45 N15-N23
      2  847 1355 (N08E19 L=163)   235     48        1 ( 817/1338- 875/1380)  33.7 E20-E17 N05-N09
      3 1726 1188 (N25W39 L=222)   547    119        1 (1687/1142-1774/1234) -18.4 W35-W44 N22-N26
      4 1825  932 (N14W55 L=237)   464     66        2 (1802/ 904-1844/ 961)  63.4 W52-W56 N11-N15
      5  860  901 (S17E06 L=177)   579    157        0 ( 782/ 851- 930/ 944)  52.1 E09-E02 S22-S13
      6  219  874 (S35E53 L=129)   596     64        0 ( 195/ 845- 237/ 911)  76.0 E55-E51 S36-S32
      7 1593  735 (S05W42 L=225)  1471    378        4 (1555/ 599-1622/ 909) -60.6 W36-W47 S12-N04
      8 1155  669 (S22W18 L=201)   365    100        2 (1108/ 644-1205/ 699)   8.1 W14-W21 S22-S20
      9  651  572 (S43E11 L=172)  1375    367        2 ( 490/ 549- 821/ 600)  18.4 E25-W02 S48-S38
     10 1492  382 (S26W54 L=236)  1225    291        0 (1366/ 355-1599/ 415)  32.0 W45-W62 S31-S20
  

    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