Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1602 1459 (N28W48 L=301)  1608    251        0 (1537/1389-1681/1555) -67.4 W43-W54 N22-N34
      2  352 1418 (N32E54 L=200)   687    102        1 ( 299/1392- 396/1446)  -8.1 E61-E47 N29-N33
      3  289 1339 (N26E58 L=196)  1536    202        1 ( 238/1307- 372/1407) -11.9 E67-E48 N24-N30
      4 1212 1339 (N23W14 L=268)   205     64        1 (1179/1321-1243/1357) -36.9 W12-W16 N21-N24
      5  681 1248 (N20E22 L=232)   416    132        1 ( 623/1205- 733/1288)  35.5 E25-E18 N17-N22
      6 1657  881 (S09W42 L=296)   812    156        1 (1629/ 829-1690/ 952)  66.8 W39-W45 S11-S04
      7 1749  859 (S11W51 L=305)   515     41        1 (1723/ 839-1777/ 886)   0.0 W48-W53 S12-S09
      8  574  564 (S23E35 L=219)   686    142        0 ( 552/ 500- 601/ 629)  71.6 E36-E32 S27-S18
      9  385  549 (S24E53 L=201)   527     81        1 ( 350/ 525- 414/ 580)  12.5 E57-E48 S25-S21
  

    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