Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1451 1611 (N25W47 L= 48)  1154    214        1 (1352/1566-1542/1639) -39.3 W41-W52 N20-N29
      2 1163 1474 (N26W22 L= 22)   552    155        1 (1105/1446-1247/1510) -40.6 W18-W25 N22-N29
      3 1377 1427 (N17W33 L= 33)   462     99        1 (1328/1405-1431/1447) -45.0 W29-W35 N14-N19
      4 1577 1380 (N09W44 L= 45)   774    176        0 (1529/1326-1608/1441) -90.0 W40-W46 N05-N13
      5  976 1281 (N20W04 L=  5)  1150    413        2 ( 820/1219-1158/1335) -39.1 E03-W13 N12-N27
      6 1372 1287 (N10W27 L= 28)   297    112        1 (1331/1249-1407/1327) -68.2 W25-W28 N07-N13
      7 1749 1223 (S04W51 L= 52)   568     88        1 (1699/1201-1787/1250) -12.5 W46-W55 S05-S02
      8 1343 1212 (N07W23 L= 24)   242     75        0 (1325/1183-1367/1243)  56.3 W21-W24 N04-N08
      9 1456  974 (S09W23 L= 24)   531    161        1 (1385/ 946-1526/1016)  -5.7 W18-W28 S10-S07
     10 1146  887 (S06W03 L=  4)   244     63        0 (1116/ 861-1170/ 908)  18.4 W01-W04 S06-S04
     11 1720  651 (S38W39 L= 40)  2622    295        7 (1654/ 528-1769/ 780)  59.0 W32-W44 S45-S29
     12  732  603 (S11E28 L=333)   631    192        0 ( 694/ 523- 792/ 656)  33.7 E33-E23 S15-S08
     13 1260  519 (S31W00 L=  1)  1634    305        4 (1209/ 378-1296/ 644)  55.5 E07-W04 S39-S23
     14 1558  543 (S39W22 L= 23)   521     59        1 (1524/ 521-1599/ 570) -23.2 W19-W26 S40-S37
     15 1424  480 (S39W10 L= 11)   559     83        3 (1378/ 463-1472/ 514) -20.6 W06-W14 S41-S36
     16  975  267 (S40E30 L=331)   831     99        1 ( 941/ 239-1028/ 305) -81.9 E32-E27 S43-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