Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1557 1558 (N19W51 L=248)   556    107        1 (1519/1525-1603/1576) -26.6 W48-W55 N18-N21
      2  797 1396 (N34E03 L=193)   264     85        0 ( 769/1359- 817/1438)  26.6 E06-W00 N31-N35
      3  541 1306 (N35E24 L=172)   631    188        3 ( 468/1257- 633/1352)   8.1 E29-E15 N33-N36
      4  165 1284 (N42E60 L=136)   572    147        1 ( 134/1265- 198/1302)   9.5 E63-E56 N40-N42
      5  435 1125 (N27E35 L=162)   309     59        1 ( 406/1103- 459/1148)  11.3 E37-E32 N24-N27
      6  383  957 (N18E42 L=154)   247     59        0 ( 359/ 932- 410/ 986) -63.4 E43-E39 N16-N20
      7 1042  532 (S22E13 L=183)   236     71        0 (1011/ 508-1076/ 560) -63.4 E13-E11 S24-S19
      8  500  448 (S14E51 L=145)  1672    432        6 ( 413/ 362- 559/ 536) -74.7 E56-E42 S20-S08
  

    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