Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  323 1158 (S10E48 L=263)   938    164        3 ( 245/1101- 381/1231)  45.0 E54-E44 S13-S04
      2  679 1030 (S12E21 L=291)   261     64        1 ( 651/1011- 713/1058) -14.0 E22-E18 S13-S11
      3  350  918 (S24E43 L=268)  1145    428        0 ( 221/ 827- 489/1051)  37.3 E54-E33 S30-S14
      4  578  991 (S16E27 L=285)   523    180        1 ( 500/ 951- 663/1046)  -9.5 E32-E20 S17-S13
      5  143  953 (S23E64 L=247)  1888    301        4 ( 116/ 860- 177/1035)  73.3 E66-E60 S28-S18
      6 1845  931 (N06W58 L=  9)  1190    196        6 (1807/ 861-1889/ 984) -22.6 W53-W65 N03-N09
      7 1035  738 (S23W06 L=318)   536    258        0 ( 968/ 644-1105/ 830)  72.9 W04-W08 S29-S16
      8 1571  664 (S16W43 L=354)   276     73        0 (1542/ 639-1599/ 684)  53.1 W40-W44 S17-S13
      9  498  636 (S40E30 L=281)   290     62        0 ( 469/ 618- 538/ 651)  15.9 E33-E26 S41-S37
     10  308  577 (S47E50 L=261)   446     47        0 ( 293/ 552- 325/ 607)  53.1 E52-E48 S48-S44
     11 1632  514 (S22W54 L=  5)   467    128        1 (1580/ 492-1678/ 555)  50.2 W50-W56 S24-S18
     12 1591  347 (S32W61 L= 12)  2350    399        4 (1464/ 313-1680/ 392)  21.3 W49-W69 S36-S27
  

    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