Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1215 1901 (N57W54 L=349)   842     70        0 (1194/1861-1234/1934)  19.3 W45-W65 N52-N59
      2 1243 1700 (N39W32 L=327)  1470    243        1 (1125/1660-1339/1734) -21.8 W22-W39 N34-N41
      3 1555 1617 (N27W52 L=348)   946    169        1 (1482/1562-1631/1658) -48.8 W47-W56 N22-N30
      4 1533 1523 (N22W45 L=340)  1205    243        1 (1473/1463-1631/1565) -48.4 W39-W50 N16-N25
      5  942 1468 (N29W03 L=298)  1633    523        1 ( 699/1439-1179/1492) -14.5 E13-W18 N24-N34
      6 1459 1457 (N19W36 L=332)   403     24        0 (1436/1430-1482/1484) -90.0 W34-W37 N17-N20
      7 1386 1401 (N17W30 L=325)   589    190        0 (1295/1381-1469/1419) -14.0 W23-W35 N14-N18
      8  603 1379 (N30E22 L=274)   334     81        0 ( 554/1358- 641/1399) -15.9 E25-E18 N27-N31
      9  369 1340 (N31E42 L=254)  1167    237        1 ( 271/1295- 459/1391) -16.7 E53-E33 N27-N35
     10  960 1250 (N15W00 L=295)  1004    370        0 ( 824/1152-1103/1330) -49.1 E07-W06 N07-N22
     11 1734 1157 (S03W49 L=344)   512     66        0 (1697/1136-1765/1180)   0.0 W45-W51 S04-S02
     12 1862 1067 (S11W61 L=356)   783    134        0 (1844/1000-1875/1141)  69.4 W57-W63 S14-S06
     13  766  967 (N02E16 L=279)   287    105        0 ( 720/ 932- 802/ 991)  18.4 E19-E13 N00-N03
     14 1604  931 (S14W35 L=330)   228     93        0 (1562/ 906-1643/ 961)   8.1 W31-W38 S15-S13
     15 1500  752 (S23W26 L=321)   387    191        0 (1417/ 709-1581/ 798)  13.0 W19-W32 S24-S21
     16 1216  758 (S18W07 L=302)   192     11        0 (1196/ 738-1233/ 774)   0.0 W05-W08 S18-S16
     17  995  722 (S16E07 L=288)   358    119        0 ( 952/ 687-1058/ 765) -54.5 E08-E03 S19-S12
     18 1050  497 (S31E08 L=287)  1205    239        1 (1012/ 391-1102/ 618)  57.3 E13-E02 S38-S24
  

    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