Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1314 1678 (N39W28 L=  7)  5211    794       10 (1013/1483-1616/1816) -30.3 W05-W46 N23-N52
      2  743 1584 (N36E17 L=322)  1871    502        2 ( 528/1549- 951/1616) -11.3 E36-E01 N32-N39
      3 1134 1575 (N32W12 L=350)   222     62        0 (1097/1560-1166/1589)   0.0 W08-W13 N31-N33
      4  448 1390 (N24E38 L=301)  2863    552        5 ( 231/1239- 596/1497) -17.7 E69-E25 N14-N34
      5 1651 1429 (N20W47 L= 26)   705     93        3 (1616/1397-1680/1475)  45.0 W43-W50 N17-N22
      6 1202 1421 (N21W14 L=353)   348    139        0 (1153/1393-1259/1469) -33.7 W11-W17 N19-N24
      7 1052 1298 (N14W03 L=342)   896    429        1 ( 906/1216-1234/1415) -35.8 E05-W13 N08-N22
      8  680 1331 (N19E20 L=319)   236    128        0 ( 632/1294- 723/1372)  39.8 E23-E17 N16-N21
      9  170 1285 (N20E67 L=272)   599     61        1 ( 158/1254- 185/1322) -76.0 E70-E63 N18-N22
     10 1667 1257 (N09W44 L= 23)   342     59        0 (1643/1234-1695/1278) -26.6 W42-W46 N07-N10
     11 1387  896 (S12W22 L=  1)  1114    440        1 (1290/ 697-1484/1060)  58.2 W15-W28 S23-S02
     12  781  853 (S11E16 L=323)   938    273        0 ( 710/ 741- 856/ 953)  45.0 E21-E10 S17-S05
     13 1590  840 (S16W36 L= 15)   569    232        2 (1522/ 748-1667/ 913)  42.3 W31-W42 S21-S11
     14 1134  719 (S21W05 L=344)   201     62        0 (1111/ 694-1164/ 741) -45.0 W03-W06 S22-S19
     15  333  541 (S27E59 L=280)   663    136        0 ( 319/ 479- 349/ 606)  45.0 E63-E54 S30-S22
     16  913  487 (S35E12 L=327)   266     77        0 ( 879/ 467- 950/ 505)  21.8 E14-E09 S36-S33
     17  583  387 (S39E43 L=296)  1023    166        0 ( 513/ 365- 669/ 426) -17.1 E48-E34 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