Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1284 1475 (N28W28 L=184)   603    123        2 (1237/1439-1341/1507)   0.0 W23-W32 N25-N29
      2  491 1267 (N31E32 L=124)  1599    459        2 ( 376/1119- 656/1408) -45.0 E42-E20 N20-N42
      3  225 1148 (N28E56 L=100)  2022    412        7 ( 180/1039- 294/1253) -60.3 E60-E49 N20-N34
      4 1633 1094 (S03W40 L=195)  1637    301        4 (1518/ 978-1748/1149)  -5.7 W29-W49 S07-N00
      5  683 1075 (N16E20 L=136)   399    117        1 ( 641/1055- 751/1096)  -8.1 E22-E15 N15-N18
      6  220 1002 (N19E58 L= 98)   501    117        0 ( 206/ 949- 245/1061) -74.1 E59-E55 N15-N22
      7 1684 1026 (S08W42 L=198)   608    175        1 (1626/ 990-1743/1049) -24.0 W37-W47 S10-S05
      8  872  909 (N03E11 L=145)  1109    217        2 ( 805/ 813- 946/1010) -81.3 E13-E07 S03-N10
      9 1704  837 (S20W41 L=196)   683    178        4 (1639/ 795-1751/ 899)  15.3 W34-W45 S22-S17
     10 1559  840 (S16W29 L=185)   474     92        1 (1544/ 784-1576/ 887)  81.9 W27-W31 S19-S12
     11  862  577 (S17E19 L=136)   909    175        3 ( 785/ 503- 907/ 626) -66.0 E23-E16 S21-S12
     12  704  283 (S35E43 L=113)  1199    223        4 ( 660/ 235- 750/ 324)  21.0 E50-E37 S38-S31
  

    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