Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1515 1745 (N38W58 L=237)   581     59        2 (1477/1725-1549/1764)   0.0 W53-W64 N36-N38
      2 1560 1488 (N20W44 L=222)   483    137        0 (1504/1443-1625/1525) -45.0 W41-W48 N16-N22
      3  464 1461 (N36E34 L=144)  1050    231        1 ( 385/1400- 572/1502)   6.7 E41-E24 N32-N38
      4 1424  984 (S08W23 L=202)   497    160        0 (1369/ 923-1491/1022)  20.6 W19-W27 S10-S05
      5 1003  684 (S19E07 L=172)   665    272        0 ( 907/ 609-1082/ 789)  31.6 E13-W00 S23-S14
      6  240  612 (S13E67 L=111)  2351    174        7 ( 197/ 530- 273/ 688)  74.7 E77-E60 S18-S07
      7  804  618 (S21E20 L=158)   290     79        0 ( 762/ 603- 847/ 632) -21.8 E22-E17 S21-S19
      8  740  583 (S22E25 L=153)   242     66        1 ( 728/ 544- 756/ 620)  59.0 E26-E23 S24-S19
      9 1256  573 (S31W08 L=186)   267     64        0 (1223/ 552-1280/ 603)  21.8 W05-W10 S31-S28
     10  531  460 (S26E45 L=133)   395     89        0 ( 511/ 423- 557/ 503)  71.6 E47-E43 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