Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  808 1477 (N26E14 L=359)   429    201        0 ( 745/1402- 842/1527)  58.0 E17-E11 N20-N29
      2  199 1323 (N18E64 L=309)  3069    220        8 ( 153/1233- 237/1408)  59.0 E72-E58 N13-N24
      3  264 1313 (N17E55 L=317)   795    168        4 ( 233/1253- 290/1379) -76.0 E57-E52 N13-N21
      4  539 1096 (N03E30 L=342)   407     98        1 ( 508/1055- 580/1140)  51.3 E32-E27 N00-N05
      5  381  764 (S17E45 L=327)  1259    326        1 ( 302/ 632- 461/ 895) -64.8 E49-E40 S25-S08
      6  157  729 (S17E73 L=299)   982     92        1 ( 117/ 693- 195/ 756) -10.0 E84-E67 S19-S15
      7 1695  616 (S28W51 L= 63)   794    131        1 (1653/ 557-1728/ 678)  74.1 W48-W53 S32-S24
      8  383  474 (S35E57 L=315)  1217    133        2 ( 321/ 436- 451/ 505)   7.4 E70-E47 S37-S33
      9 1282  343 (S49W22 L= 34)  1503    227        4 (1175/ 309-1383/ 388) -10.3 W11-W33 S51-S44
     10 1572  371 (S46W52 L= 65)   469     38        0 (1547/ 355-1600/ 386)   9.5 W48-W56 S47-S44
     11 1050  338 (S48W00 L= 13)   406     88        0 (1019/ 318-1095/ 363)   0.0 E02-W04 S50-S46
     12  660  226 (S57E48 L=324)   546     64        1 ( 631/ 211- 697/ 245)   4.8 E53-E41 S58-S55
  

    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