Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1735 (N51W39 L=282)   570     65        1 (1310/1701-1358/1766)  63.4 W36-W42 N47-N53
      2 1848 1319 (N16W70 L=313)  1528    142        1 (1786/1277-1900/1364) -23.6 W63-W80 N12-N19
      3  474 1157 (N14E37 L=206)   323     61        0 ( 442/1136- 497/1178) -21.8 E39-E34 N13-N15
      4  678  875 (S03E22 L=220)   624    256        1 ( 580/ 827- 797/ 939) -23.2 E28-E14 S06-N00
      5  913  699 (S15E09 L=234)   225     58        0 ( 892/ 669- 937/ 722)  53.1 E10-E06 S17-S13
      6 1565  624 (S25W37 L=280)  1099    235        3 (1490/ 549-1626/ 691)  42.0 W31-W41 S29-S20
      7  304  656 (S17E57 L=186)   569     63        1 ( 281/ 630- 340/ 683)  14.0 E60-E52 S19-S15
      8  291  594 (S22E62 L=181)   491     55        1 ( 277/ 562- 315/ 621)  33.7 E65-E58 S23-S19
      9  719  530 (S26E24 L=219)   400    141        0 ( 677/ 479- 783/ 570)  31.0 E28-E18 S29-S23
  

    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