Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1181 1857 (N59W47 L= 56)   957     99        2 (1138/1822-1223/1889)   7.4 W37-W60 N56-N61
      2  997 1808 (N60W21 L= 29)   329     41        0 ( 971/1784-1019/1829)  12.5 W16-W25 N58-N62
      3 1525 1645 (N34W55 L= 64)   609    158        1 (1456/1605-1598/1666) -35.0 W49-W59 N30-N37
      4 1489 1543 (N29W44 L= 53)   717    157        1 (1408/1511-1554/1570) -33.7 W39-W48 N26-N32
      5  350 1278 (N30E47 L=322)  2502    361        3 ( 242/1170- 467/1414) -29.5 E66-E35 N22-N39
      6  593 1323 (N31E25 L=344)   705    192        1 ( 576/1225- 614/1405)  84.8 E26-E23 N24-N35
      7 1638 1341 (N14W47 L= 56)   441    113        0 (1602/1305-1675/1366)  15.9 W43-W50 N12-N15
      8 1215 1238 (N16W15 L= 24)   819    331        2 (1064/1202-1348/1302) -34.5 W06-W22 N11-N22
      9 1747 1145 (N01W51 L= 59)   347     68        1 (1718/1122-1785/1166)   0.0 W47-W54 N00-N01
     10  462  948 (N09E37 L=332)   475     83        2 ( 445/ 899- 482/1001)  71.6 E38-E35 N06-N12
     11 1833  920 (S15W56 L= 65)  1656    196        4 (1764/ 860-1892/ 992)  14.0 W48-W64 S17-S12
     12  735  686 (S10E23 L=346)   680    138        1 ( 710/ 621- 761/ 748)  60.3 E25-E20 S13-S05
     13 1036  634 (S17E05 L=  3)   716    249        0 ( 942/ 593-1107/ 710)  20.0 E11-W00 S19-S13
     14  904  687 (S12E12 L=356)   220     69        0 ( 866/ 673- 938/ 702) -36.9 E14-E10 S13-S10
     15  240  647 (S08E66 L=303)   724     74        0 ( 225/ 608- 260/ 683)  45.0 E69-E62 S10-S05
     16  355  547 (S15E57 L=311)   978    247        2 ( 290/ 499- 404/ 608) -76.0 E62-E50 S18-S10
     17  294  538 (S16E65 L=303)   716     64        2 ( 270/ 508- 320/ 571) -53.1 E68-E62 S17-S13
     18  563  358 (S30E48 L=320)  1508    199        5 ( 504/ 294- 608/ 433)  24.4 E62-E39 S35-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