Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1455 1724 (N33W59 L=137)  1128     90        2 (1406/1686-1498/1758) -63.4 W56-W63 N29-N35
      2 1823 1106 (S12W55 L=133)   439     46        0 (1804/1084-1847/1134) -59.0 W53-W57 S14-S09
      3 1734  969 (S16W43 L=121)  1272    218        1 (1710/ 869-1752/1088)  72.9 W39-W46 S22-S09
      4  581 1015 (N17E27 L= 51)   224     40        1 ( 553/ 993- 606/1033) -33.7 E28-E25 N15-N18
      5  331  857 (N13E50 L= 28)   543     82        2 ( 306/ 810- 349/ 891) -68.2 E51-E47 N10-N15
      6  185  722 (N07E72 L=  6)   620     60        2 ( 158/ 699- 208/ 745) -36.9 E76-E67 N05-N08
      7 1528  405 (S48W18 L= 96)   367     75        0 (1508/ 376-1555/ 437)  38.7 W16-W21 S50-S46
      8  885  421 (S25E26 L= 51)   279     91        0 ( 833/ 413- 931/ 432) -26.6 E29-E23 S26-S23
      9  880  272 (S36E36 L= 42)   960    199        0 ( 807/ 222- 932/ 349)  18.4 E45-E27 S38-S31
     10 1404  304 (S52W04 L= 82)   427     52        0 (1376/ 290-1440/ 318) -15.9 W01-W08 S54-S50
  

    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