Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1699 (N21E53 L=280)   361     62        1 ( 551/1664- 587/1728) -33.7 E56-E50 N18-N22
      2 1721 1500 (N46W39 L= 12)  1472    185        2 (1657/1426-1778/1567) -32.5 W30-W45 N42-N49
      3 1648 1360 (N32W31 L=  4)   868    140        2 (1600/1296-1679/1420) -29.1 W25-W34 N28-N35
      4 1581 1229 (N21W27 L=  1)   295     64        0 (1556/1198-1597/1258) -33.7 W25-W28 N19-N22
      5 1898 1064 (N24W62 L= 35)   556     72        0 (1883/1029-1913/1097) -31.0 W59-W65 N21-N25
      6  943  939 (S14E02 L=331)   199     79        0 ( 909/ 917- 969/ 966)  63.4 E03-W00 S15-S11
      7  783  831 (S24E09 L=325)   294    155        0 ( 709/ 805- 855/ 850)  36.9 E12-E04 S27-S21
      8 1253  579 (S24W28 L=  1)   768    255        0 (1146/ 539-1370/ 612)  42.3 W22-W34 S29-S19
  

    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