Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  537 1520 (N25E36 L=226)   544    143        0 ( 508/1461- 568/1581)  82.9 E37-E34 N21-N29
      2 1684 1428 (N25W45 L=307)  1917    323        0 (1570/1291-1768/1537) -51.3 W36-W51 N16-N31
      3  417 1476 (N22E45 L=217)   585    141        0 ( 360/1450- 468/1494) -15.3 E50-E39 N20-N23
      4 1363 1421 (N22W20 L=282)   212     42        0 (1340/1403-1392/1436)  26.6 W18-W22 N20-N22
      5 1615 1258 (N13W38 L=299)  1312    212        1 (1540/1188-1712/1338) -31.6 W32-W45 N09-N17
      6  729 1295 (N11E19 L=243)   199     40        0 ( 711/1272- 747/1316) -56.3 E20-E18 N09-N12
      7 1163 1259 (N10W07 L=269)   486    157        0 (1089/1236-1250/1289)  18.4 W03-W12 N08-N12
      8 1813 1156 (N09W55 L=316)   349     56        1 (1795/1123-1828/1185)  76.0 W52-W56 N06-N10
      9 1849 1114 (N06W59 L=320)   914    194        4 (1808/1074-1900/1154)  10.3 W54-W65 N03-N09
     10  152 1056 (S03E65 L=196)  1094    183        1 ( 127/ 969- 176/1140)  90.0 E69-E62 S08-N01
     11  370 1012 (S06E43 L=219)   921    188        1 ( 340/ 925- 433/1075)  60.9 E45-E37 S11-S02
     12  531  964 (S09E31 L=231)   325    107        1 ( 495/ 925- 582/ 997)  33.7 E33-E27 S11-S07
     13  719  798 (S19E19 L=243)  1076    294        0 ( 621/ 691- 860/ 870)  34.5 E25-E09 S25-S14
     14 1707  789 (S14W48 L=310)   509    105        1 (1663/ 759-1740/ 831)   0.0 W44-W51 S15-S11
     15 1535  608 (S27W38 L=300)  1211    218        2 (1446/ 554-1600/ 711)  54.0 W31-W41 S30-S19
     16  590  545 (S35E30 L=231)   665    212        1 ( 536/ 474- 653/ 604)  39.3 E36-E24 S40-S31
     17 1753  561 (S28W59 L=321)   995    143        1 (1696/ 510-1799/ 591)  40.6 W54-W63 S31-S25
  

    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