Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1581 1655 (N46W54 L=255)  1540    243        3 (1493/1613-1691/1692)  -6.3 W43-W71 N42-N49
      2 1044 1556 (N37E00 L=201)   420     87        0 (1000/1537-1088/1575)   9.5 E03-W03 N35-N38
      3  336 1501 (N30E60 L=141)   709     63        1 ( 311/1468- 372/1528)  23.2 E63-E56 N27-N31
      4 1109 1486 (N33W05 L=206)   334     64        1 (1081/1445-1132/1516) -68.2 W02-W06 N29-N34
      5 1106 1342 (N23W04 L=206)   397    115        0 (1084/1299-1156/1374)  45.0 W03-W07 N20-N24
      6  189 1173 (N08E65 L=137)   766    133        0 ( 176/1111- 207/1241) -53.1 E67-E61 N04-N12
      7  780 1193 (N12E16 L=185)   221     64        0 ( 747/1168- 815/1219) -21.8 E18-E13 N11-N13
      8 1506 1098 (N08W31 L=232)  1133    297        2 (1420/ 961-1566/1199) -60.3 W25-W35 N00-N14
      9 1662 1058 (N06W43 L=244)   349     53        1 (1631/1037-1694/1084) -31.0 W40-W45 N04-N07
     10 1739  817 (S09W51 L=252)  1171    163        2 (1698/ 739-1776/ 897)  73.3 W48-W54 S13-S03
     11  212  840 (S12E61 L=140)   494     74        1 ( 187/ 821- 243/ 860)  20.6 E64-E56 S13-S10
     12  196  730 (S20E67 L=134)   855     74        0 ( 168/ 692- 218/ 766)  33.7 E74-E62 S22-S16
     13 1646  591 (S24W48 L=249)  1510    180        2 (1584/ 517-1704/ 653)  52.1 W42-W52 S28-S19
     14  710  601 (S25E20 L=181)  1043    227        2 ( 636/ 566- 805/ 650)  -4.8 E25-E13 S27-S21
     15  503  593 (S27E37 L=165)   367     98        0 ( 455/ 577- 541/ 618)  14.0 E40-E32 S27-S24
     16 1192  218 (S55W22 L=223)  2410    248        5 (1080/ 187-1315/ 250)   9.1 W10-W36 S59-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