Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  948 1794 (N52E24 L=223)   510    104        0 ( 902/1777- 990/1814)   0.0 E29-E20 N50-N53
      2  308 1511 (N20E63 L=184)  1571    230        3 ( 277/1458- 345/1575) -53.1 E67-E57 N16-N24
      3  404 1383 (N13E48 L=199)   296    126        0 ( 366/1346- 435/1406) -20.6 E52-E44 N11-N14
      4  190 1222 (N01E67 L=180)   566     56        0 ( 173/1194- 205/1250)  90.0 E69-E64 N00-N02
      5 1231 1182 (N13W11 L=258)   498    209        0 (1136/1166-1327/1221)   0.0 W04-W16 N11-N14
      6 1342  950 (N00W20 L=268)  1138    262        1 (1306/ 860-1414/1066) -70.0 W17-W26 S04-N07
      7  309  862 (S20E48 L=199)   505    140        0 ( 271/ 823- 364/ 912) -26.6 E51-E42 S21-S16
      8  400  870 (S18E40 L=207)   330     73        0 ( 379/ 833- 427/ 902)  68.2 E41-E37 S20-S15
      9  809  879 (S11E11 L=236)   419    129        0 ( 743/ 864- 882/ 891)  12.5 E15-E06 S13-S09
     10 1127  861 (S08W08 L=256)   219     76        0 (1090/ 838-1170/ 877)  11.3 W06-W11 S09-S07
     11 1345  611 (S21W27 L=274)  1944    561        6 (1162/ 460-1486/ 782)  62.5 W18-W33 S33-S08
     12 1047  742 (S17W05 L=253)   589    303        1 ( 942/ 716-1206/ 762)  18.4 E01-W15 S19-S14
  

    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