Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  456 1632 (N20E53 L=275)  1157    270        1 ( 426/1579- 507/1699) -26.6 E60-E46 N16-N24
      2  397 1404 (N04E49 L=280)   399     76        1 ( 374/1367- 414/1442) -45.0 E50-E46 N02-N06
      3 1039 1279 (N08E04 L=324)   603    219        1 ( 975/1218-1127/1311)  14.0 E07-W01 N03-N09
      4  396 1072 (S15E41 L=288)   246     80        0 ( 361/1044- 426/1095)  38.7 E43-E38 S17-S13
      5 1329  852 (S10W21 L=350)   200     57        0 (1291/ 832-1361/ 875) -11.3 W18-W23 S11-S09
      6 1115  828 (S16W09 L=338)   218     60        0 (1089/ 801-1142/ 867)  68.2 W08-W10 S18-S13
      7  880  767 (S25E03 L=325)   647    312        2 ( 738/ 701-1026/ 815)  34.5 E11-W05 S31-S19
      8  861  205 (S66W18 L=346)   621    117        2 ( 837/ 172- 888/ 233) -22.6 W12-W24 S69-S63
  

    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