Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1407 1739 (N35W53 L= 33)   345     77        0 (1372/1716-1442/1759) -78.7 W51-W54 N32-N37
      2 1592 1650 (N24W61 L= 41)  1151    142        2 (1510/1629-1658/1670) -24.8 W54-W68 N21-N27
      3 1223 1566 (N30W28 L=  7)   488    138        0 (1160/1524-1270/1587) -39.8 W24-W30 N26-N32
      4 1467 1566 (N23W44 L= 24)   504    159        0 (1392/1557-1548/1576) -24.4 W39-W50 N20-N25
      5 1452 1507 (N20W40 L= 20)   382    101        0 (1410/1478-1508/1532) -38.7 W37-W42 N17-N22
      6 1787 1102 (S11W52 L= 31)   712    134        1 (1724/1080-1850/1128)  -8.7 W45-W58 S13-S09
      7  473  981 (N13E35 L=305)   235     39        0 ( 457/ 954- 489/1010)  45.0 E36-E33 N11-N14
      8 1563  559 (S39W26 L=  5)  1148    285        1 (1472/ 465-1671/ 667)  24.0 W17-W36 S42-S33
      9  459  344 (S25E64 L=275)  1480    134        3 ( 412/ 300- 515/ 391) -26.6 E74-E57 S27-S21
     10  691  321 (S31E45 L=295)   939    173        2 ( 611/ 282- 763/ 363) -33.7 E50-E40 S34-S26
     11  553  291 (S30E60 L=279)  1570    156        1 ( 530/ 219- 584/ 355)  23.5 E76-E52 S35-S25
     12  778  213 (S42E48 L=291)   541     59        0 ( 764/ 178- 796/ 252)  31.0 E53-E43 S44-S38
  

    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