Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  662 1679 (N45E35 L=128)  2799    474        7 ( 540/1599- 807/1731)   4.8 E44-E20 N38-N50
      2 1292 1630 (N44W21 L=184)   408     77        0 (1276/1591-1313/1660)  59.0 W18-W23 N41-N46
      3 1000 1615 (N42E04 L=159)   356     64        1 ( 976/1586-1023/1648)  68.2 E06-E02 N39-N44
      4 1418 1542 (N37W29 L=193)   308     55        1 (1393/1523-1447/1568)  11.3 W26-W31 N36-N39
      5 1000 1397 (N26E03 L=161)   272    105        1 ( 943/1358-1050/1427) -26.6 E06-W00 N23-N28
      6  511 1246 (N15E35 L=128)   735    205        3 ( 451/1188- 555/1315) -45.0 E40-E31 N11-N18
      7  608 1247 (N15E28 L=135)   255     98        1 ( 563/1232- 654/1257)   0.0 E31-E24 N14-N15
      8  646 1190 (N12E25 L=139)   333     57        1 ( 627/1159- 671/1224)  90.0 E26-E23 N09-N13
      9  522 1152 (N09E33 L=130)   766    255        2 ( 455/1104- 636/1215) -24.8 E38-E25 N06-N12
     10 1608  573 (S25W44 L=208)   451    117        0 (1570/ 544-1642/ 599) -20.6 W40-W48 S27-S23
     11  321  486 (S36E58 L=106)   724    227        1 ( 291/ 458- 345/ 514)  29.1 E63-E54 S39-S33
     12 1432  418 (S37W35 L=198)  1206    207        1 (1363/ 395-1518/ 462)  -4.4 W28-W41 S38-S33
     13 1050  397 (S39W04 L=168)   798    225        2 ( 944/ 376-1164/ 420)  10.6 E04-W13 S41-S36
  

    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