Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1142 1584 (N40W05 L=340)   431     69        2 (1096/1566-1179/1604)   0.0 W01-W08 N38-N41
      2  668 1475 (N28E29 L=306)   619    207        0 ( 591/1430- 775/1512) -14.0 E36-E20 N26-N30
      3 1651 1392 (N29W46 L= 21)  1200    245        0 (1579/1304-1721/1476) -50.7 W40-W51 N23-N34
      4 1434 1425 (N30W27 L=  2)   309     76        1 (1402/1403-1470/1449) -11.3 W24-W29 N28-N31
      5  490 1312 (N16E39 L=296)  1488    320        3 ( 392/1215- 592/1421) -30.1 E50-E30 N10-N22
      6  964 1304 (N19E06 L=329)   623    209        2 ( 871/1251-1037/1351) -22.6 E12-W00 N16-N21
      7  787  724 (S18E14 L=321)   331     76        1 ( 745/ 706- 824/ 746)   0.0 E16-E11 S19-S16
      8 1445  694 (S16W30 L=  5)   385     76        1 (1417/ 660-1474/ 726)  68.2 W28-W31 S18-S13
      9 1703  654 (S18W53 L= 28)   408     48        0 (1683/ 628-1727/ 676)  56.3 W51-W55 S19-S16
     10  294  542 (S35E64 L=271)   736     92        1 ( 249/ 523- 336/ 567)   9.9 E78-E55 S37-S32
  

    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