Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1381 1496 (N28W18 L=318)   316     79        0 (1347/1454-1405/1520) -45.0 W15-W20 N25-N29
      2  631 1388 (N13E29 L=271)  1145    244        1 ( 579/1301- 695/1479) -50.7 E33-E23 N07-N18
      3 1805 1156 (N13W52 L=352)  2102    424        3 (1719/1049-1914/1217) -10.8 W44-W65 N07-N16
      4 1446 1109 (N04W25 L=325)   320    102        2 (1395/1068-1490/1138)  -9.5 W21-W27 N02-N05
      5 1339 1028 (S02W19 L=319)   286    104        0 (1300/ 997-1388/1068) -15.9 W15-W22 S03-N00
      6  152  962 (S15E64 L=236)  1081    210        3 ( 109/ 918- 199/1003)  -9.5 E70-E58 S17-S12
      7  675  891 (S17E20 L=280)   229     71        1 ( 633/ 877- 713/ 907)  18.4 E23-E17 S17-S15
      8  183  872 (S20E60 L=240)   937    124        1 ( 144/ 840- 213/ 898) -23.2 E64-E56 S22-S18
      9 1721  776 (S11W50 L=350)   434    104        0 (1691/ 739-1757/ 808) -23.2 W47-W54 S12-S09
     10 1357  614 (S26W27 L=327)  1779    396        2 (1181/ 539-1556/ 746)  39.0 W16-W37 S32-S15
     11 1041  510 (S37W08 L=308)   447    119        2 ( 972/ 489-1094/ 533)  20.6 W03-W11 S39-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