Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  402 1657 (N24E70 L=230)   807     57        1 ( 379/1638- 428/1686)  -6.7 E80-E63 N22-N25
      2  593 1642 (N27E47 L=253)   482    123        2 ( 560/1605- 632/1671) -78.7 E50-E43 N23-N28
      3 1203 1446 (N27W02 L=302)   940    327        2 (1077/1345-1321/1506) -12.5 E07-W11 N23-N30
      4 1383 1424 (N30W15 L=315)   275    104        0 (1351/1379-1417/1470) -39.8 W11-W17 N27-N32
      5 1877 1188 (N27W65 L=  5)  2372    226        1 (1843/1082-1901/1290)  85.9 W59-W69 N20-N34
      6  445 1063 (S11E37 L=262)   536     99        0 ( 395/1036- 482/1090)  33.7 E41-E34 S13-S09
      7 1765  990 (N12W51 L=351)   764    105        1 (1738/ 935-1790/1036)  74.1 W48-W53 N08-N15
      8 1448  559 (S20W39 L=339)  1986    352        3 (1365/ 436-1579/ 654)  80.0 W35-W45 S28-S11
      9  471  554 (S44E33 L=267)   340     51        1 ( 457/ 526- 491/ 580)  68.2 E34-E30 S46-S41
     10  919  532 (S34W05 L=305)   636    221        1 ( 825/ 502-1023/ 559)  30.3 E01-W11 S37-S29
  

    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