Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  981 1803 (N41E33 L=337)   408     50        0 ( 950/1788-1018/1815)  45.0 E35-E30 N38-N42
      2  605 1659 (N19E48 L=322)  1340    312        2 ( 580/1523- 638/1797) -45.0 E59-E40 N10-N29
      3 1073 1759 (N40E24 L=346)   515     67        1 (1044/1732-1103/1785) -15.9 E27-E20 N38-N42
      4  709 1606 (N19E38 L=332)   795    185        1 ( 683/1519- 730/1689) -39.8 E44-E32 N13-N23
      5  475 1618 (N14E56 L=314)   388     67        0 ( 456/1584- 494/1648) -78.7 E58-E54 N11-N16
      6 1524 1446 (N34W19 L= 29)   962    144        2 (1495/1377-1553/1510) -63.4 W15-W22 N29-N38
      7  792 1385 (N08E23 L=347)  1214    304        3 ( 689/1288- 913/1472)  65.0 E26-E18 N00-N15
      8  253 1236 (S13E58 L=312)  1036    221        1 ( 167/1210- 353/1257)  14.0 E69-E48 S15-S10
      9 1015 1166 (N01E04 L=  6)   965    271        3 ( 922/1127-1097/1224)  16.7 E10-W00 S01-N06
     10 1624 1052 (N13W35 L= 45)   328    138        0 (1588/ 987-1651/1113) -35.5 W30-W38 N10-N15
     11  596 1065 (S15E26 L=344)   309     79        1 ( 557/1041- 629/1089)  45.0 E28-E23 S17-S13
     12 1509 1013 (N07W28 L= 38)   467    101        0 (1491/ 946-1524/1061) -63.4 W25-W30 N04-N10
     13  224  870 (S36E57 L=313)  1515    224        2 ( 140/ 793- 298/ 944)  23.2 E76-E47 S41-S29
     14 1185  837 (S12W14 L= 24)   193     36        0 (1166/ 815-1209/ 853)  71.6 W13-W15 S13-S10
     15  942  712 (S26W05 L= 15)  1344    347        1 ( 790/ 661-1100/ 776)  47.0 E03-W12 S32-S17
     16 1514  628 (S13W42 L= 52)   349    134        0 (1459/ 611-1581/ 644)  18.4 W38-W47 S14-S11
     17 1363  496 (S24W39 L= 49)   336     98        0 (1310/ 486-1414/ 505)  29.7 W35-W42 S26-S22
     18 1237  337 (S38W42 L= 52)  1652    283        2 (1129/ 271-1358/ 414)  79.4 W38-W46 S45-S29
  

    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