Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1017 1873 (N61E31 L=111)  1387    125        5 ( 965/1831-1072/1915)   0.0 E45-E21 N56-N65
      2  350 1613 (N24E72 L= 69)  1188     87        2 ( 321/1571- 374/1652) -56.3 E77-E68 N21-N27
      3  798 1544 (N28E26 L=116)   475    233        1 ( 746/1505- 877/1611)   0.0 E30-E19 N25-N32
      4  726 1413 (N18E26 L=115)   228     73        0 ( 706/1385- 754/1448) -38.7 E28-E23 N16-N20
      5  723 1260 (N09E22 L=119)   280    121        0 ( 674/1241- 776/1290)   0.0 E26-E18 N08-N09
      6 1722 1185 (N21W46 L=188)   775    138        3 (1694/1107-1742/1252) -77.5 W43-W48 N16-N25
      7  379  805 (S25E42 L=100)   881    228        1 ( 329/ 716- 450/ 870)  47.5 E47-E35 S31-S19
      8  526  736 (S27E29 L=113)   485     94        2 ( 499/ 705- 562/ 778) -45.0 E31-E26 S28-S24
      9  816  690 (S24E08 L=134)   249     99        1 ( 782/ 650- 842/ 726)  63.4 E09-E06 S27-S21
     10 1260  562 (S25W24 L=166)  2533    601       12 (1026/ 450-1479/ 720)  48.7 W11-W34 S36-S11
     11 1712  664 (S10W56 L=197)   365     41        1 (1684/ 642-1736/ 683)  14.0 W52-W58 S11-S08
  

    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