Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  730 1444 (N16E28 L= 90)   417    134        2 ( 684/1400- 768/1486) -33.7 E32-E25 N13-N18
      2  950 1390 (N18E13 L=105)   710    157        2 ( 875/1359-1023/1419)  29.1 E17-E08 N14-N20
      3 1618 1069 (N14W36 L=155)   319    105        0 (1605/1020-1637/1124) -59.0 W34-W38 N11-N16
      4  718  814 (S21E14 L=104)   370     88        1 ( 690/ 774- 755/ 843)  78.7 E16-E12 S23-S18
      5  623  732 (S28E19 L= 99)   604    149        1 ( 556/ 705- 689/ 761)  31.0 E24-E14 S31-S25
      6  387  705 (S36E40 L= 78)   359     34        0 ( 371/ 674- 404/ 730) -71.6 E41-E37 S37-S33
      7 1487  456 (S25W48 L=166)   786    108        1 (1441/ 425-1536/ 491)  20.6 W45-W53 S27-S22
      8  570  451 (S50E22 L= 96)   373     66        0 ( 529/ 437- 606/ 461)  20.6 E26-E18 S51-S48
      9  804  392 (S48W01 L=120)   842    140        2 ( 727/ 373- 878/ 414)  14.0 E05-W07 S50-S45
     10 1383  356 (S34W47 L=165)   408     89        0 (1347/ 327-1425/ 390)  80.5 W45-W49 S36-S30
     11 1174  170 (S55W54 L=172)   523     62        0 (1160/ 136-1188/ 209) -18.4 W46-W64 S57-S51
  

    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