Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  655 1701 (N59E23 L= 23)   468     25        0 ( 631/1675- 676/1723)  76.0 E26-E19 N57-N61
      2  520 1651 (N57E40 L=  6)  2259    217        5 ( 429/1591- 614/1702)   6.6 E55-E27 N52-N61
      3  850 1570 (N44E02 L= 43)   646    152        3 ( 783/1547- 926/1596) -31.0 E07-W03 N41-N47
      4 1071 1532 (N37W14 L= 59)   395    142        1 (1001/1515-1137/1548)  -5.2 W08-W19 N35-N38
      5 1184 1319 (N21W15 L= 61)   673    194        0 (1118/1276-1266/1343) -18.4 W10-W20 N18-N23
      6  773 1100 (N15E14 L= 31)   581    185        0 ( 731/1021- 799/1181) -84.3 E16-E12 N10-N20
      7  286 1058 (N18E52 L=353)   839    181        1 ( 262/ 969- 313/1126) -78.7 E54-E49 N12-N22
      8 1820  861 (S19W55 L=101)  1146    220        0 (1772/ 778-1867/ 970)  45.0 W50-W61 S23-S13
      9 1579  710 (S23W32 L= 78)  1167    246        1 (1482/ 651-1641/ 799)  24.8 W24-W37 S25-S18
     10 1133  679 (S16W01 L= 47)   593    153        1 (1105/ 610-1175/ 737)  63.4 E01-W04 S20-S12
     11 1217  525 (S28W04 L= 50)   544    112        2 (1178/ 480-1251/ 571) -66.8 W01-W05 S32-S24
  

    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