Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  536 1713 (N45E51 L=192)   567    105        1 ( 507/1689- 569/1737)  -8.1 E55-E47 N42-N46
      2 1598 1586 (N43W46 L=288)   615    150        1 (1554/1561-1641/1606)  16.7 W40-W50 N40-N44
      3 1313 1359 (N25W18 L=260)   254     78        0 (1284/1337-1356/1373)  21.8 W15-W20 N23-N25
      4 1608 1176 (N14W39 L=281)   359     68        1 (1581/1153-1635/1212) -45.0 W36-W40 N13-N16
      5 1792 1085 (N09W55 L=298)   851    161        1 (1764/1007-1814/1162) -73.3 W52-W57 N04-N14
      6 1871 1008 (N05W65 L=308)   709    137        1 (1853/ 963-1900/1064) -49.4 W62-W70 N01-N08
      7  991  887 (S06E01 L=241)   841    204        1 ( 946/ 777-1014/ 983) -72.9 E04-W00 S13-N00
      8 1407  598 (S23W29 L=271)   254     74        0 (1377/ 583-1443/ 623) -21.8 W26-W31 S23-S21
      9  393  536 (S34E47 L=196)   891    131        0 ( 356/ 514- 443/ 570)  -9.5 E50-E41 S35-S31
     10  288  522 (S37E59 L=183)   746    191        1 ( 254/ 508- 319/ 544)  12.5 E64-E55 S38-S34
  

    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