Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1599 1683 (N28W70 L= 21)  1468    143        4 (1542/1625-1653/1725) -90.0 W64-W74 N23-N31
      2 1375 1526 (N27W38 L=348)   321    108        0 (1333/1481-1415/1562) -74.1 W36-W38 N23-N30
      3 1631 1513 (N18W57 L=  7)  1016    174        1 (1554/1468-1703/1547) -45.0 W51-W61 N13-N22
      4 1506 1518 (N22W47 L=357)   269     44        0 (1494/1483-1518/1544)  36.9 W44-W48 N20-N23
      5  847 1455 (N37E00 L=310)   278     50        0 ( 828/1424- 867/1477)  76.0 E01-W01 N34-N38
      6 1026 1254 (N20W06 L=316)   192     49        0 (1015/1219-1039/1284)  63.4 W04-W07 N18-N22
      7 1761 1163 (S04W52 L=  3)   428     79        0 (1739/1120-1788/1195)  23.2 W48-W55 S06-S02
      8  199 1075 (N26E62 L=248)   525     52        1 ( 183/1046- 218/1110)  76.0 E64-E59 N23-N27
      9  274 1036 (N22E53 L=258)  1735    369        1 ( 169/ 939- 418/1107)   6.1 E67-E39 N17-N26
     10 1572  990 (S08W32 L=343)   349     48        0 (1553/ 958-1589/1022)  63.4 W30-W33 S10-S06
     11 1846  851 (S26W56 L=  7)  1793    183        3 (1787/ 784-1915/ 913)  14.0 W48-W68 S28-S22
     12 1718  870 (S20W42 L=352)   718    190        1 (1624/ 839-1771/ 911) -12.1 W33-W47 S21-S16
     13 1089  834 (S05E01 L=310)   236     49        0 (1067/ 813-1110/ 857) -71.6 E01-W00 S06-S03
     14 1520  750 (S21W24 L=335)   336    106        0 (1475/ 718-1568/ 773)   0.0 W20-W28 S22-S20
     15  507  701 (S01E40 L=271)   334     92        0 ( 496/ 654- 520/ 747)  45.0 E41-E37 S03-N01
     16 1149  548 (S24E04 L=307)   793    214        1 (1086/ 453-1212/ 641)  45.0 E10-W01 S29-S19
     17 1403  459 (S38W11 L=322)   364     67        0 (1384/ 422-1417/ 497)  71.6 W09-W12 S41-S35
     18 1004  250 (S43E27 L=284)  2755    290        7 ( 898/ 137-1065/ 365)  26.6 E49-E15 S52-S34
  

    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