Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1607 1566 (N43W28 L=152)   557     91        1 (1573/1532-1636/1602) -45.0 W24-W31 N40-N45
      2  317 1364 (S02E56 L= 68)   385     76        0 ( 294/1333- 346/1396)  45.0 E56-E53 S03-N00
      3  177 1114 (S18E65 L= 59)  3999    395       11 (  97/ 939- 246/1274)  49.4 E79-E58 S29-S08
      4 1639 1217 (N20W33 L=157)   268     78        0 (1622/1175-1655/1258) -59.0 W31-W35 N17-N22
      5 1535 1179 (N14W26 L=150)   571    158        0 (1503/1105-1556/1252) -58.0 W22-W28 N10-N18
      6 1875  862 (N08W64 L=188)  1963    158        4 (1806/ 803-1929/ 920)  31.0 W56-W72 N03-N12
      7 1396  722 (S14W30 L=153)   636    191        0 (1333/ 651-1466/ 792)  74.7 W27-W32 S19-S08
      8 1332  740 (S15W25 L=149)   250     99        0 (1293/ 708-1370/ 775)  68.2 W23-W26 S17-S12
      9  869  686 (S30E01 L=123)   830    218        1 ( 797/ 646- 966/ 734)  21.8 E06-W04 S33-S24
     10 1187  440 (S36W30 L=153)   406    112        0 (1135/ 403-1231/ 464)  -6.3 W25-W34 S37-S35
  

    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