Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1039 1775 (N51W34 L= 67)   324     52        1 (1018/1750-1058/1795)  45.0 W31-W37 N49-N53
      2 1186 1765 (N45W43 L= 76)   840     93        2 (1142/1729-1234/1794) -45.0 W39-W47 N42-N47
      3 1699 1539 (N14W65 L= 98)  1052    119        3 (1642/1523-1763/1559) -18.4 W58-W74 N11-N16
      4 1132 1374 (N23W16 L= 49)   364     71        0 (1101/1342-1155/1409)  21.8 W13-W18 N21-N25
      5  883 1209 (N21E03 L= 30)   951    356        2 ( 808/1097- 951/1372)  38.7 E10-W05 N16-N28
      6 1788 1310 (S00W60 L= 93)   551    116        1 (1757/1252-1806/1361)  68.2 W57-W61 S03-N02
      7 1219 1150 (N09W14 L= 47)   272     66        0 (1196/1125-1237/1182)  45.0 W12-W15 N07-N10
      8  288 1041 (N25E50 L=343)  1592    335        6 ( 257/ 909- 328/1152) -86.2 E53-E47 N16-N31
      9  489 1077 (N23E32 L=  1)   601    170        1 ( 405/1032- 554/1104) -31.0 E38-E27 N18-N25
     10 1671 1069 (S09W40 L= 73)   317     58        2 (1651/1034-1697/1097)  76.0 W38-W42 S11-S07
     11 1877  944 (S24W59 L= 92)   656    109        2 (1844/ 901-1894/ 997)  29.7 W54-W62 S25-S21
     12 1555  862 (S17W27 L= 60)  1161    340        2 (1409/ 836-1726/ 885) -23.5 W17-W40 S22-S12
     13  712  526 (S14E33 L=  0)   866    265        0 ( 663/ 421- 755/ 630)  38.2 E40-E26 S19-S08
     14  614  435 (S18E44 L=349)   324     55        0 ( 586/ 421- 648/ 450) -14.0 E45-E41 S19-S16
     15  587  336 (S24E53 L=340)   521    120        0 ( 536/ 317- 616/ 371)   0.0 E56-E48 S25-S21
     16  497  321 (S23E64 L=329)  1195    108        5 ( 467/ 288- 537/ 365)  50.2 E70-E59 S25-S19
     17  790  316 (S30E39 L=354)   720    135        0 ( 757/ 258- 838/ 355)  24.4 E44-E33 S33-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